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Quantitative behaviors of shock-induced dislocation nucleation are investigated by means of molec- 
ular dynamics simulations on fee Lennard- Jones solids: a model Argon. In perfect crystals, it is 
found that the Hugoniot elastic limit (HEL) is a linearly decreasing function of temperature: from 
near-zero to melting temperatures. In a defective crystal with a void, dislocations are found to 
nucleate on the void surface. Also HEL drastically decreases to 15 percent of the perfect crystal 
when the void radius is 3.4 nanometer. The decrease of HEL becomes larger as the void radius 
increases, but HEL becomes insensitive to temperature. 

PACS numbers: 62.20.Fe,62.50.+p 



Mechanical properties of shock-loaded solids are im- 
portant to materials science since they can reveal relia- 
bility of materials under extreme conditions. Shocks in 
solids generate high pressure almost instantaneously, and 
create extremely high deformation (strain) rates. Yield- 
ing phenomena realized at high strain rates are usually 
quite different from those of quasistatic ones. Notably, 
Rohde has found that the HEL is almost independent of 
temperature ; Kanel and his coworkers have found that 
the HELs of Al and Cu single crystals are increasing func- 
tions of temperature '^l . These experiments are in appar- 
ent contrast with quasistatic deformations, where yield 
strength considerably decreases with increasing temper- 
ature. However, contrary to Rohde and Kanel, stain- 
less steel shows that the HEL is a decreasing function 
of temperature These puzzling results provide us 
with intriguing problems. The microscopic foundations 
of plastic deformation are dominated by the dynamics of 
dislocations. Unfortunately, the nature of dislocation dy- 
namics is very complicated and yet to be fully understood 
from theoretical point of view. Therefore, it is reasonable 
to decompose the problem into three essential ingredients 
of dislocation dynamics: nucleation, mobility, and mul- 
tiplication. As the first stage of this line of thought, the 
properties of dislocation nucleation are investigated in 
this paper. 

In a perfect crystal, dislocations nucleate with the help 
of thermal fluctuations. Holian and Lomdahl have stud- 
ied the emergence of stacking faults initiated by partial 
dislocation emission in a shocked perfect crystal 3 . Re- 
cently, Tanguy et al. |^ have performed extensive simu- 
lations to confirm that there exists a critical size of dis- 
location loop below which nucleated loops are energeti- 
cally unstable and eventually annihilate; this is remines- 
cent of droplet nucleation in a metastable gas or liquid. 
Although these two studies are enlightening, the role of 
temperature, which is known to be important to disloca- 
tion nucleation, has not been considered. In some stud- 
ies, temperature plays no role in activating dislocation 
nucleation ^J. However, dislocation nucleation seems to 
be anyhow a fluctuation-assisted process in the sense that 
no nucleation is observed in a molecular dynamics simu- 



lation of a perfect crystal at zero temperature j^. It is 
important to understand temperature effects on disloca- 
tion nucleation to clarify the controversial experimental 
results I, 2, 3]. 

This paper discusses temperature effects on dislocation 
nucleation through molecular dynamics simulations. It 
turns out that the HEL is a linearly decreasing function of 
temperature from zero to melting temperatures. Then a 
defect is introduced into the system as a preexisting void, 
at which dislocation nucleation is found to be enhanced. 

The present simulations use a simple Lennard- Jones 
(LJ) potential U{r) = 4e [(cr/r)^^ _ (^a/rf] with a cutoff 
length of 4.0(7. To avoid a discontinutity at the cutoff, the 
potential is suitably shifted. Note that simulations with 
shorter cutoffs such as 2.6a resulted in a visible decrease 
of HEL, as much as 90 percent. Hereafter, we let a = 3.41 
Aand e = 1.65 x 10^^ [J] in order to make a quantita- 
tive comparison with experiments on solid Argon. The 
Lennard-Jones potential has been well-tested in terms of 
mechanical and thermal properties of Argon such as the 
elastic constants and the melting temperature. Hence it 
is reasonable to adopt the LJ potential as the flrst com- 
putational attempt to understand the anomalous HEL 
behavior of fee metals 0- 

In our system the axes of x, y, and z are taken to coin- 
cide with (100), (010), and (001) directions, respectively. 
Shocks propagate along the [001] direction. As usual, 
periodic boundary conditions are applied to x and y di- 
rections. The whole system consists of 100 x 100 x 200 
unit cells. Shock waves in the present simulations are 
generated by pushing a piston of inflnite mass into the 
still target. The velocity of the piston is denoted by Up. 
We call the piston velocity Up throughout this paper. 

In this work, the HEL is defined as the longitudinal 
stress (azz) above which dislocation nucleation is ob- 
served. Note that this definition of HEL differs from the 
one that is based entirely on Shockwave behavior: the 
shock front decomposes into an elastic precursor and a 
plastic wave. However, these two definitions are identi- 
cal because the decomposition is also observed in molecu- 
lar dynamics simulations once dislocations are emitted in 
shocked solids. Longitudinal stress is calculated from the 
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FIG. 1: Temperature dependence of Hugoniot elastic limits of 
a perfect crystal (circles) and of a defective crystal containing 
a void of 1.71 nm radius (squares). The dashed line represents 
experimental data for the shear modulus [^. The melting 
temperature is approximately 83 K. 



Hugoniot relation; a^z = PoUsUp, where po and Us denote 
initial density and shock velocity, respectively. Note that 
the initial pressure is zero. The longitudinal stress (Tzz is 
approximately proportional to the shear stress which 
is directly responsible for dislocation emission. 

Figure. ^ shows the nearly linear dependence of HEL 
on temperature. The HEL at 77 K, close to the melt- 
ing temperature (83 K) 7] , is approximately half that 6 
K. This is clear evidence that dislocation nucleation is 
a thermally assisted phenomenon. Also, this tempera- 
ture dependence is reminiscent of elastic constants. The 
dashed line in FIG. ^ represents the shear modulus of 
solid Argon I3| . We can see that the temperature depen- 
dences of the HEL and of the shear modulus are qual- 
itatively the same for perfect crystals. That is, HEL is 
approximately half the shear modulus regardless of tem- 
perature. (Indeed, the ratio of HEL to shear modulus is 
slightly greater than 0.5 for T < 24 K.) The critical strain 
also decreases linearly from 0.138 (at 6 K) down to 0.118 
(at 77 K) as the temperature increases. Numerical data 
are shown in Table Note that the temperature depen- 
dence of the shear modulus of Argon has been calculated 
by molecular dynamics simulation using L J potential . 

The temperature dependence of HEL obtained by the 
present simulation is opposite to the experiment on fee 
metals 0. But this is not a contradiction, for the ex- 
periment involves real crystals where various defects and 
impurities preexist. Plastic deformations of real mate- 
rials are mainly governed by the motion of preexisting 
dislocations, while in a perfect crystal only the disloca- 
tion nucleation is considered here. Hence, the mobility of 
preexisting dislocations must be analyzed for further un- 
derstanding of temperature dependence of shock-induced 
plasticity. This subject will be discussed elsewhere. 

Strikingly different results are found when we turn to 
defective crystals, concentrating on the effects of a void. 



TABLE L Critical piston velocity, Hugoniot elastic limit, and 
critical strain for the perfect crystal as functions of initial 
temperature. 
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FIG. 2: Four partial dislocation loops nucleate on the void 
surface. Void radius is 2.4 nm and piston velocity is 70 m/s. 
Atoms are colored according to the number of the nearest 
neighbors; green, red, and blue atoms have 10, 11, and 13, 
resp ectively. Atoms of 12 nearest neighbors are not presented 
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A void is introduced in a model crystal as a blanked 
spherical region. The radius is varied from 0.34 nm to 
3.41 nm. With a void, dislocations nucleate exclusively 
on the void surface as shown in FIGS. [21 El and^ 

The Hugoniot elastic limit decreases drastically in the 
presence of a void. Square symbols in FIG. ^ represent 
the HEL for a defective crystal with a void of 1.71 nm 
radius. With this size of void, the HEL is approximately 
20 percent of the perfect crystal's. Note that the HEL of 
a defective crystal also linearly decreases with increasing 
temperature up to 42 K. 

When r > 42 K, approximately half the melting tem- 
perature, no dislocation is nucleated and a void collapses 
after the passage of shock front. Although this phe- 
nomenon is similar to "hotspot" formation this is 
not the case since the piston velocity used here is not 
strong enough to cause a jet which is essential to hotspot. 
The collapse of the void observed in the present situation 
rather postulates that voids are unstable to perturbations 
above half the melting temperature • 

Below 42 K, the HEL and the critical piston velocity 
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FIG. 3; Thirteen picoseconds after FIG. El Partial dislocation 
loops are extendend and Lomer-Cottrel sessile dislocation is 
formed. 




FIG. 4: Sideview of FIG. El Two slip planes on the void 
surface are activated. 



for the defective crystal with a void become insensitive 
to temperature. For the case of r = 1.02 nm, the critical 
piston velocity is 103 ± 5 m/s at 6 K, dropping only to 
95 ± 5 m/s at 42 K, while for the perfect crystal the 
corresponding values are 289 ± 1.6 m/s and 213 ± 1.6 
m/s, respectively. The absolute decreases are 8 m/s for 
the defective crystal and 76 m/s for the perfect crystal. 
Then the relative changes are 8/103 ~ 0.08 and 76/284 ~ 
0.26, respectively. The difference between the perfect 
and the defective crystals is remarkable in terms of both 
absolute and relative decreases. Thus defects dominate 
thermal fluctuations in the enhancement of dislocation 
nucleations. In experiment, the microstructure of the 
specimen is important. 

The void-size dependence of HEL is shown in Fig. 
|S1 At small r, the HEL decreases as followed by 
a crossover to r^°-^ at Tc = 1.3 ± 0.3 nm. Note that the 
transition radius rj, is approximately the width of a dou- 
ble kink and is also very close to the critical loop radius 
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FIG. 5: Log-log plot of HEL as a function of void radius. 
Solid and dashed lines are proportional to r~°'^ and r~^, re- 
spectively. Initial temperature is 6 K. 

TABLE II: Comparison of simulation and experimental shock 
velocities as functions of piston velocity. (Simulation data 
involves a perfect crystal). Initial temperature and density 
for the simulation are 64 K and 1.63 g/cc, respectively, while 
they are 75 K and 1.65 g/cc in the experiment. 
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1.5 nm above which nucleated dislocation loops grow sta- 
bly The quantitative explanation of this intriguing 
transition is not clear at this point. 

To compare the present simulation with experiment, 
the Hugoniot values {up versus Ug) of Argon for both this 
simulation and experiment |l3j | are presented in TABLE 
Hi! We see that the simulation values exceed the experi- 
mental ones by approximately 15 percent. This overesti- 
mation is comparable to that of Bclonoshko's simulations 
using a Buckingham potential [l^- His deviation from 
the experimental value exceeds 10 percent, especially for 
small piston velocities up to 1.2 km/sec. Dubrovinsky 
has attributed this overestimation to the polycrystalline 
nature of the experimental specimen [l5(. Since both 
Belonoshko's and the present simulations deal only with 
crystals of (100) orientation, the deviation from the ex- 
perimental result is not unreasonable. 

Finally, two quantitative uncertainties in the present 
simulation are remarked. One is underestimation of 
stacking fault energy by the use of LJ potential. It 
is known that stacking fault energy vanishes when one 
adopts short-range, two-body intermolecular potentials 
[l5l |. In this regard, it is possible that the numerical val- 
ues of the HEL obtained by the present simulations are 
underestimates. But, at least, qualitative results such as 
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the temperature dependence or the void effect are not 
so influenced by the nature of two-body potentials. In 
addition, a large cut-off, 4. Oct, was used to reduce this 
effect. 

The other uncertainty invloves periodic boudary con- 
ditions. The system size (57 x 57 x 114 nm) is small 
enough that a dislocation loop may interact with itself 
through the periodic boudary conditions. It is impossible 
to remove this uncertainty from molecular dynamics sim- 
ulations at this point. However, the strain around a nu- 
cleated dislocation loop propagates at the sound velocity 
(about 900 m/s for transverse mode), taking 45 picosec- 
onds to cross the simulation cell. Therefore, this efect 
may not seriously influence the present results, which in- 
volve exclusively nucleation processes that occur within 
10 picoseconds. 

To summarize, molecular dynamics simulations on 
both perfect and defective fee crystals are performed from 



near-zero to melting temperatures. It is found that the 
Hugoniot elastic limit is about half the shear modulus 
and decreases linealy with increasing temperature. The 
critical piston velocity and critical strain also decrease 
in the same manner. When a void is introduced in the 
crystal, dislocations are found to nucleate on the void 
surface. The critical piston velocity decreases drastically 
as the void radius r increases. The decrease is nonlinear 
and shows a crossover from to r"*^'^ around r ~ 1.3 
nm, which is approximately the width of a double kink. 
The Hugoniot clastic limit becomes insensitive to tem- 
perature in the presence of a void. That is, the effect of 
nucleation at void dominates that of thermal fluctuations 
in shock-induced dislocation nucleation. 

The author gratefully acknowledges helpful discussions 
with H. Kaburaki, F. Shimizu, Y. Satoh, N. Nita, and H. 
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